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Abstract 

The existence of compact dimensions which are accessible only to gravity represents 
an intriguing possible solution to the hierarchy problem. At present the strongest 
constraint on the existence of such compact Gravity-Only Dimensions (GODs) 
comes from SN 1987a. Here we report on the first self-consistent simulations of 
the early, neutrino-emitting phase of a proto-neutron star which include energy 
losses due to the coupling of the Kaluza-Klein modes of the graviton which arise 
in a world with GODs. We compare the neutrino signals from these simulations to 
that from SN 1987a and use a rigorous probabilistic analysis to derive improved 
bounds for the radii of such GODs. We find that the possibility that there are two 
compact extra dimensions with radii larger than 0.66 fim is excluded at the 95% 
confidence level — as is the possibility that there are three compact extra dimensions 
larger than 0.8 nm. 
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1 Introduction 



The possibility that there exist large, compact, space-like dimensions accessi- 
ble only to gravity, but not to the standard model particles, has recently been 
suggested as a solution to the hierarchy problem [1,2]. If such compact spatial 
dimensions do exist then the fundamental Planck scale could be close to the 
electro-weak scale, with gravity becoming comparable in strength to other in- 
teractions at energies of order a few TeV. At present, the strongest constraint 
on the existence of compact Gravity-Only Dimensions (GODs) such as these is 
derived from the SN 1987a neutrino signal. The argument, put simply, is that 
if too much of the available energy is radiated to the GODs then the neutrino 
signal seen on terrestrial detectors would have differed markedly from that ac- 
tually observed. This argument allowed Arkani-Hamed et al. (ADD) to place 
order-of-magnitude bounds on the radii of the GODs [2]. Subsequently, an 
attempt to improve the analysis of Ref. [2] was made by Cullen and Perelstein 
in Ref. [3], who found R < 0.3 /im {R < 0.4 nm) for 2 (3) compact extra di- 
mensions. Further improvements to the calculation of the crucial Kaluza-Klein 
(KK) graviton emissivities were made in Ref. [4]. These calculations tended 
to confirm that the original estimate of Ref. [2] was indeed within an order of 
magnitude of the correct number, although they resulted in bounds a factor 
of two or three weaker: R < 0.71 fim (R < 0.85 nm). However, as indicated 
in Ref. [3], to obtain truly accurate bounds on the size of the GODs — and, by 
inference, the fundamental Planck scale — one would have to incorporate the 
contribution to the emissivities for the graviton emission to the extra dimen- 
sions into a numerical code for PNS evolution. In contrast, Refs. [3,4] both 
used the simple criterion suggested by Raffelt [5], which requires only that the 
emissivity (energy radiated per unit mass and per unit time) of some exotic 
particle be less than 10^^ ergs/g/s. 

This bound is based on the time-scale of several seconds over which the 
handful of electron-type anti-neutrino events from SN 1987a were detected 
in Kamiokande [6] and 1MB [7]. These detections confirmed the standard sce- 
nario of core-collapse Supernovae: neutron stars are born in the aftermath of 
successful supernova explosions, as the stellar remnant becomes gravitation- 
ally decoupled from the expanding ejecta [8-12]. The ambient conditions in 
a newly born neutron star, also called a proto-neutron star (PNS), are so 
extreme — densities of order lO^"' g/cm"^ and temperatures between 10 and 50 
MeV — that the neutrino mean-free path is much smaller than the stellar ra- 
dius. Gonsequently, the enormous gravitational binding energy gained during 
the stellar collapse (~ 2 — 3 x 10^^ ergs) is stored inside the PNS and slowly 
released by neutrino diffusion on a time scale of tens of seconds. This pic- 
ture meshes nicely with the long duration of SN 1987a's neutrino signal, since 
the neutrinos must have diffused, and not free-streamed, out of the PNS. In 
this context, any weakly-interacting particle which couples to nucleons and 
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can freely stream out of the star easily competes with neutrinos as a means 
of transporting energy away from the stellar interior. Such an additional heat 
sink would accelerate the cooling and remove a significant fraction of the bind- 
ing energy. The late-time neutrino signal, which is fueled by the heat stored in 
the core, is particularly sensitive to these effects, and it would be drastically 
altered if the couplings to exotic particles were strong enough. 

It is this physics that facilitates the use of the supernova neutrino signal 
to place bounds on the amount of energy lost as exotic radiation. Stringent 
bounds have been placed on the properties of a number of novel, weakly- 
interacting, light-particle species. Initially, these inchidcd limits on the neu- 
trino magnetic moment, the mass of the axion, the strength of right-handed 
neutrino interactions and the sterile- neutrino production rates [13-18]. Subse- 
quently, several authors have employed Raffelt's bound on the total emissivity 
due to these exotic energy- loss mechanisms — or analogous bounds on the lumi- 
nosity [19] — as a simple means of enforcing the constraint provided by the SN 
1987a neutrino signal. However, since the emission rates are strong functions 
of temperature, the bounds on exotic-particle couplings which are deduced 
from such arguments will depend crucially on the fiducial value of temper- 
ature at which the exotic-particle emissivity is calculated. Not surprisingly, 
since the temperature of the newly-born neutron star varies with both posi- 
tion and time, the fiducial temperature is rather loosely defined, with values 
between 10 and 70 MeV appearing in the literature. Furthermore, exotic cool- 
ing mechanisms will lower the temperature inside the PNS core, and so the 
appropriate fiducial temperature is really a function of the emissivity of these 
processes. The strong temperature dependence of these emissivitics ultimately 
means that these ambiguities can easily change the resulting bound by more 
than an order of magnitude. 

Thus in order to establish more accurate bounds it is necessary to perform 
detailed simulations of the effect of exotic radiation on the neutrino signal. 
So far such simulations have only been carried out for the case of axion ra- 
diation [18,20]. These simulations tend to confirm the picture outlined in the 
previous two paragraphs. Therefore in this work we report on the first self- 
consistent simulations of the PNS core which test the effect of GODs on the 
SN 1987a signal. In §2 we present a summary of the early stages of the life 
of a neutron star and the effect of additional energy loses. In order to make 
these simulations as realistic as possible we use the most recent calculations 
of the KK-graviton emissivity [4] , which differ by roughly a factor of five from 
earlier, more schematic, calculations [3]. 

Another important drawback of former analyses is the lack of a well-established 
selection criteria to disregard a given model. As emphasized in Ref. [21], the 
sparseness of the SN 1987a data calls for a Likelihood analysis. In §3 we out- 
line the basics of the likelihood formalism which allows for marginalization 
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over uncertain parameters, as well as the derivation of a bound with a clear 
probabilistic interpretation and unambiguously-defined confidence levels. Our 
results, and a discussion of them, are presented in §4. 



2 PNS Evolution with Additional Sources of Energy Loss and SN 
1987a signal 



Core-collapse supernovae are incredibly rich physical systems, which include 
many complicated phenomena. Indeed there is still much debate about basic 
features like the details of the explosion mechanism [22,23]. What is, how- 
ever, now well-established is that the birth and subsequent evolution of the 
PNS located in the inner region of the supernova provides the main source 
of neutrinos during the first several seconds after core bounce. Nearly all of 
the binding energy gained during the collapse is trapped inside the PNS and 
this energy is radiated by neutrino diffusion over several tens of seconds [8- 
11]. These time scales and the associated neutrino luminosities arc affected 
by several pieces of physics, including the total mass of the PNS, and both 
the nuclear equation of state (EOS) and neutrino opacities at supra-nuclear 
density. None of these quantities can be determined in a model-independent 
approach, but earlier studies have explored the model sensitivities by using 
different EOSs and the associated self-consistently-calculated opacities [11]. 
These studies indicate that the critical physical parameter that determines 
neutrino luminosity in the first several seconds is the PNS total mass. The 
uncertainties associated with the properties of dense hadronic matter do not 
affect these luminosities greatly, provided that the opacities are calculated 
consistently with the EOS, because important feedbacks between the EOS 
and the opacities tend to reduce the differences. The protoncutron star evolu- 
tion code employed for this study has been described in detail in Refs. [11,12]. 
Details regarding the neutrino opacities employed can be found in Ref . [24] . In 
a recent paper the rates for the emissivities of KK-gravitons in dense matter 
were computed in a model independent way using low-energy theorems which 
relate the emissivities to the well measured nuclcon-nucleon cross sections [4]. 
In a world with n = 2 or n = 3 compact GODs the KK-graviton emissivity 
of neutron-star matter can be fitted to within 5% accuracy over a range of 
temperatures from 15 to 30 MeV and a range of densities from 0.5no to 3no, 
where uq = 0.16 fm"'^ is the nuclear saturation density. The result, which 
includes contributions from nn, pp, and np collisions is: 

dE fnn\ / T 

dr = «"(^)(l0l*v) ^(^"■^'' MeV/baryon/s (1) 
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with p2 = 5.42 , 02 = 5.1 x lO^Rj, ps = 6.5, 03 = 1.4 x lO^^Rl Here, ub is 
the baryon number density, T is the temperature, and the size of the extra 
dimensions, i?„, is given in mm [^. The emissivity depends on the composition 
of matter through the function 

x(X„,X,)=X^ + Xp2 + 4X„X,(^^^^j , (2) 

where is the neutron fraction, and Xp = (1— X„) is the proton fraction. The 
weak, residual, temperature dependence in the function x arises because the 
np cross section decreases as approximately l/\/E^, where Ecm is the center 
of mass energy of the colliding nucleons, while the nn cross section is roughly 
independent of energy, for the energies of interest. Over the range of values 
of proton fraction and temperature encountered in the PNS the function x 
changes by less than 10% from the case of pure neutron matter, when x = 1- 
Thus, the total emissivity is a weak function of the proton fraction Xp, and 
in our simulations we assume, for simplicity, that X„ = 1. 

In this paper we will put limits on permitted values of the coefficients 02 and 
as, assuming the temperature-dependencies T^'^^ and T^'^. Thus, ultimately 
our results are valid for the coupling to any exotic particle which results in 
these temperature-dependencies. By placing limits on the coefficients a„, we 
are led, in this case, to bounds on i?„ — the radii of the extra dimensions — 
but the simulations, philosophy, and likelihood formalism discussed here can 
be applied in a much broader context. The count rate in the detectors is 
sensitive to the neutrino luminosities and spectrum. The procedure employed 
to obtain the count rates is described in detail in earlier work [11,12]. In 
addition, here we include the effects due to the neutron-proton mass difference 
on the anti-neutrino absorption cross sections and correct the fiducial mass of 
the 1MB detector to 6 ktons [25]. Both lead to better agreement between the 
two detectors. 

The equilibrium diffusion approximation employed in the numerical code pro- 
vides a fair description of the total neutrino luminosity and its time structure. 
However, additional assumptions are required if an anti-neutrino spectrum 
is to be derived from this information. In particular, to obtain an average 
anti-neutrino temperature — denoted here by Tp^ — we assume that the neutri- 
nos leaving the protoneutron star are thermalized inside a neutrinosphere and 
have a Fermi-Dirac spectrum with zero chemical potential. The neutrinosphere 
was chosen to be at an optical depth of 2/3 in accordance with the usual def- 

^ In the non-degenerate limit Ref. [4] derived an analytic formula with p2 = 5.5, 
P3 = 6.5, and 02 = 5.6 x lO^i?!; '^3 = 1-5 x lO^^iZ^. One might perhaps be justified 
in using these formulae, but instead we have chosen to use more accurate fits to the 
full result of Ref. [4]. 
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Fig. 1. Comparison of the results of our simulations for a PNS with a baryonic mass 
of M = 1.6 Mq to the data from SN 1987a. N(t) denotes the integrated number of 
counts. We show results from a model without graviton emission (solid line), as well 
as the results with graviton emission for 02 = 0.01 MeV/baryon/s (short-dashed 
line) and 02 = 0.1 MeV/baryon/s (long-dashed line). 

inition. However, it must be borne in mind that the neutrino mean-free paths 
are energy dependent, and thus the concept of a neutrinosphere is only ap- 
proximate. In the next section we will study the effect of those uncertainties 
related to this definition of the neutrino temperature by exploring the sensitiv- 
ity of the count rates and the deduced bounds to changes in our prescription. 
In Fig. 1 we show the number of accumulated counts as a function of time, 
N{t), in Kamioka (left panel) and 1MB (right panel). The results shown are for 
a baryonic mass of 1.6 Mq (with a final gravitational mass of 1.46 Mq), and 
an anti-neutrino temperature defined as above. They include the case without 
KK-graviton emission, and two cases with radiation to the GODs included, 
namely: 02 = 0.01 MeV/baryon/s and a2 = 0.1 MeV/baryon/s. As argued in 
the introduction, KK-gravitons steal energy from the core and thereby dra- 
matically suppress the late-time neutrino emission. Note that the early-time 
neutrino emission is not strongly influenced by the existence of GODs because 
these neutrinos are emitted from the lower density and lower temperature re- 
gions of the star where the graviton emissivity is small. However, after about 
2-4 seconds, neutrino losses provoke a strong compression and heating of the 
star [8,10,11], which activates the KK-graviton emissivity. Consequently, the 
total number of neutrino counts, which clearly depends on the total amount 
of energy radiated in neutrinos, drops as 02 is increased. 

This total number of counts also has a significant dependence on the PNS 
mass. Indeed, the mass of the PNS is generally the most important factor in the 
neutrino signal [11]. However, it is worth noting that, if a sufficiently important 
additional sink of energy is included in the PNS simulation, the time-structure 
of the signal actually becomes less sensitive to the precise value of the PNS 
mass. This happens because more massive stars attain higher temperatures, 
which in turn results in enhanced graviton emission. Consequently, much of 
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the additional binding energy reservoir in a more massive star is carried out 
by exotic particles, and the resultant neutrino luminosity is not dissimilar to 
that from PNSs of lower mass. This was already noticed in Ref. [20] for the 
axion case. 



3 From simulations to bounds: the desirability of a likelihood anal- 
ysis 

Nevertheless, this "mass-independence" only arises for large values of the cou- 
plings a2 and a^. More typically we face the problem of trying to extract 
bounds for these couplings from the neutrino data in the face of significant un- 
certainties in the PNS mass M and the anti- neutrino temperature Tp^ . There- 
fore a key question is how to perform a consistent analysis of the results and 
derive a bound on the coupling of the exotic particles. Studies similar to ours 
for the case of axion emission [18,20], based their investigations on either the 
time required to accumulate 90% of the counts (Ai(90%)) or the total number 
of counts detected N^°^. 

Not surprisingly, both A'^*"* and At (90%) decrease with increasing R2. The 
variation in the total number of counts is moderate and very sensitive to 
the variation of all three parameters studied: 02, Tp^ and M. Thus, it makes 
more sense to study an observable that is relatively insensitive to Tp^ and 
M. At(90%) is such a quantity. However, it is still difficult to move from the 
values of At(90%) to a confidence level for a bound on R2. One might argue 
that a simple criterion is to disregard those radii R2 which lead to At (90%) 
lower than one half of the value for the a„ = case. Using this criterion we 
find that: 

R2 < (0.4 ... 0.9) fxm (3) 

for PNS masses in the range M = (1.5 ... 2.0) Mq. However, there is no 
guarantee that the result for a„ = is the "correct" answer. The performance 
of a model that badly overestimates the number of counts in the a„ = 
case might well be improved by the presence of GODs. So, if we adopt an 
alternative approach and disregard those values of R2 which lead to a time 
scale that is less than one half of the value found in the SN 1987a data we 
find: 

R2 < (0.2 ... 0.4) //m . (4) 
These, however, are rather arbitrary criteria. There is really no statistical 
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basis [26] for judging what constitutes a successful reproduction of the exper- 
imental data, and thus no way to define confidence levels. 

In light of these issues we seek an analysis tool that will allow us to av- 
erage over parameters which are not well-known (such as the mass of the 
PNS and TpJ and also give a clear criterion that allows us to decide what 
values of the parameter space are excluded by the data. To this end we 
now pursue a "likelihood" analysis. The goal of this analysis is the proba- 
bility that the coupling a„ is smaller than some specified value a°, given the 
data, and a certain set of model assumptions. We denote this probability by 
prob(a„ < a^\{data}, I), where the I indicates the set of physical assumptions 
which underly the simulations discussed in the previous section. Note that the 
assumptions denoted by / do not include information about M or Tp^. Our goal 
is to relate prob(a„ < a^\{data}, I) to a function we can calculate directly from 
our model. That function is the Likelihood function CD{{data}\an, M, Tp^, /), 
which is the the probability that the data actually taken at a detector D (KII 
or 1MB) arises given a set of model assumptions (/) and certain specific values 
for the key quantities M, Tp^, and a„ 0. Assuming that the data obeys Poisson 
statistics, the likelihood for the neutrino signal in a single detector can then 
be expressed in terms of the model prediction for the count rates at the times 
where an event happened and the total number of events actually seen [27] 

CD{{data}\an, M,Tp^, I) = 



when Nf) is the total number of observed neutrino arrivals, and 
are the times at which neutrinos actually arrived in the detector D. Noit) 
is the total number of neutrinos that the model with assumptions /, PNS 
mass M, anti- neutrino temperature Tp^, and exotic-coupling a„ predicts will 
have arrived in the detector D up until the time t. At can be any interval 
small enough that the probability of detecting more than one count in any 
one bin can be taken to be negligible. Ultimately it will be absorbed into an 
overall normalization constant. Note that in order to get the full answer for 
the likelihood C the product of the likelihoods for KII and 1MB must be taken 
(i.e., C = Ckii X Cimb). 

The likelihood function provides us a with quantitative tool to compare models 
of PNS evolution. Since we treat not only the coupling to extra dimensions, 
an, but also M and Tp^, as parameters, the likelihood function is a function in 

^ We could, of course, have allowed for explicit variation of other parameters of 
the PNS simulation, and not simply included them in /, as we do when we cal- 
culate CD{{data}\an, M,Tp^, I). The formalism can easily be extended in this way 
if it becomes clear that some other physical parameter affects the neutrino signal 
strongly. 
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a three-dimensional space. This function has a minimum at a baryonic mass of 
M = 1.5AfQ, 0,2 = and Tj^^ = 1.1T°. T° denotes the "reference anti-neutrino 
temperature" defined above using the optical-depth prescription. We can now 
assess all other models by looking at the log of the ratio of their likehhood of 
a particular model to this "most hkely" model: 

C{{data}\an,M,Tp^,I) \ 
C{{data}\0, l.SM©, I.IT^", I) J ^ ' 



q{an,M,T,-J = - log 



The function q is then a function in this same three-dimensional space. 

For the case of two GODs we discuss first the situation where these dimen- 
sions have zero radius, and so 02 = 0. Contours of q in the resulting two- 
dimensional X — M plane are shown in the left panel of Fig. (2). Here the 
horizontal axis is the temperature normalized to the reference anti-neutrino 
temperature, x — Tp^/T°, which is varied to explore the sensitivity of the re- 
sults to deviations of the neutrino temperature from that obtained using the 
optical depth prescription. We see that although M = 1.5Mq, T = 1.1T° is 
indeed a minimum of the function g it is a rather shallow minimum: vary- 
ing the neutron-star mass and the anti-neutrino temperature over the range 
considered here does not reduce the hkelihood greatly. Ultimately this reflects 
the weakness of the constraint that the SN 1987a data provides for these 
parameters. 

The situation is rather different as we move away from 02 = 0. In the middle 
panel we display similar likelihood contours in the M — log(a2) plane for 
the case Tp^ = T°, while the right panel shows contours in the x — log(a2) 
plane for the baryonic mass M = 1.6Mq. (The pattern of contours shown 
is typical, and does not change significantly if different values of Tp^ and M 
are chosen.) These panels show that the hkelihood function decreases rapidly 
for a2 > 10~^ — regardless of the values of the poorly- known parameters M 
and Tp^. Such large values of 02 are essentially two orders of magnitude less 
likely than the "most likely" model. This is a contrast to smaller values of 
02, where the differences in likelihood are comparable to those seen in the 
M — Tp^ plane. Thus statements about the most hkely value of 02 in this 
smaller-a2 regime will be sensitive to the ill-constrained information on these 
neutron-star parameters. ConscqTiently we will not make any such statements 
here. However, the two rightmost panels give us confidence in our ability to 
derive a bound on 02 — as opposed to a most-likely value — since it is clear that 
certain values of this coupling can be well-excluded, completely independent 
of details of the PNS modehng. 

With the likelihood function already in hand, the derivation of such a bound is 
quite a simple matter. A natural ansatz for the probability of a„ to be smaller 
than some specified a° is 
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Fig. 2. Contour levels of the function q defined in Eq. (6): {x — M) left panel; 
(M — log(a2)) middle panel; {x — log(a2)) right panel. 
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prob(an < al\{data}J) = N J da J dM g{M) j dTpJ{%^,M) 



X £({dato}|a„,M,rp,,7). (7) 

In other words, we now integrate (or "marginalize" ) over all possible values of 
M and Tj>^, using appropriate weight functions. These functions are f{Tp^, M), 
which is the probability for the neutrinosphere to be at a temperature Tp^ 
given a PNS mass, M, and g{M), which is the probability for that PNS mass 
to occur. The overall constant N is defined by the normalization condition 
prob(a„ > 0|{ciaia},/) = 1 . 

A rigorous derivation of this expression can be given by employing Bayesian 
statistics [28]. In the course of this derivation the weight functions / and g ac- 
quire the aforementioned probabilistic interpretations. Note that the structure 
of Eq. (7) is easily generalized: for every parameter not well determined from 
other sources one must sum over the parameter space with the appropriate 
relative weight. In our case we will assume that the mass of the PNS hes within 
a certain range [M^j„, M^aa,], but that within that range all values are equally 
probable, i.e., we write g{M) oc 6{M — Mmin)d{Mmax — M). Meanwhile, the 
central value of the anti-neutrino temperature is calculated according to the 
optical-depth prescription, and T^^ is then assumed o be Gaussian distributed 
with a 10% width, independent of time. It might be argued that a potential 
systematic error exists in the choice of anti-neutrino temperature, and thus a 
different weight function should be used. To check this we also used a function 
/ which was constant over a range of temperatures that differed by ± 20% 
from the central value. This altered the bound obtained by less than 10%. 
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4 Improved bounds and discussion 

Equation (7) may now be used to find bounds on 02 and 03. Defining the bound 
at the 95% confidence level to be that value of a° which obeys prob(a„ < 
a^\{data}, I) = 0.95, we find that, for the n = 2 case 

= 0.022 MeV/baryon/s and thus R2 < 0.66 /xm. (8) 

A similar analysis for the n = 3 case differs in its details, but yields the same 
overall picture. The bound (again, at the 95% confidence level) is: 

al = 0.0077 MeV/baryon/s and thus R3 < 8 x 10-Vm. (9) 

These are within 10% of the bounds derived in Ref. [4], where the emissivity 
(1) was derived and the Raffelt criterion applied using ns = uq and a fidu- 
cial temperature chosen rather arbitrarily to be 30 MeV. The good agreement 
between this work and Ref. [4] might be thought to suggest that simple crite- 
ria based on energetics are an entirely satisfactory. However, analyses of our 
simulations suggest that the actual fiducial temperature in the 02 = case 
was 35-40 MeV, and the density was as large as 4no for the most massive 
stars. Application of the same criterion at these temperatures and densities 
could have resulted in a bound on 02 that was too stringent by a factor of ten. 
In contrast, the fiducial temperature in our simulations with 02 > 0.01 was 
generically less than 20 MeV, which leads to a bound on 02 that is too weak 
by an order of magnitude. The advantage of our analysis of the SN 1987a data 
is that it obviates any need to make a subjective judgment about the value of 
the "real" fiducial temperature. Such an analysis confirms that, for n = 2 and 
n = 3, the supernova bound on the size of extra dimensions found from the 
neutrino signal of SN 1987a is significantly more stringent than any collider 
bound which will be obtained in the foreseeable future. It is also, roughly, a 
factor of hundred more stringent than the bound derived by means of recent 
tests of sub millimeter gravity [29] for the case of two extra dimensions. We 
regard this astrophysical information as a robust bound, since our 95% confi- 
dence level is fairly insensitive to uncertainties in the neutron-star mass and 
the anti-neutrino temperature [Q These bounds are specific to scenarios where 
the extra dimensions are fiat or weakly warped [32]. They come about because 
of the existence of a large number of low-lying KK modes with energy less the 
characteristic temperature of the supernova core, and thus do not apply to 
extra-dimensional models where such low-lying KK modes are absent [33]. 

^ Arguments from cosmology [30,31] may provide tighter constraints on R2 and R3, 
but there is much we presently do not understand about cosmology in the presence 
of GODs. 
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We have not studied the extent to which changes in the neutrino mean-free 
path or the equation of state of strongly-interacting matter at high density 
influence our results. Recent work on the influence of many-body correlations 
on the neutrino-transport mean-free path indicates that the mean-field value 
for this quantity may be a factor of two or three too small [34]. However, 
as discussed in Ref. [20] in the axion-radiation case, increases in the mean- 
free path actually strengthen the bound on exotic-particle couplings, since the 
shorter neutrino diffusion times that result make it harder to reconcile the 
late-time neutrino signal with the presence of exotic radiation mechanisms. 
Consequently we view the bound derived here as one that is conservative in 
regard to its assumptions about neutrino mean-free paths. 

A potentially more significant physics uncertainty is our lack of knowledge 
about the specific heat of matter at high density. Changes in this function 
might lead to a large effect on the bounds for R2 and R3, since they will 
produce different interior temperatures for the PNS, and the KK-graviton 
emissivity is a strong function of temperature. However, any changes in the 
PNS's interior temperature will also affect the neutrino diffusion time since 
neutrino mean- free paths are strongly temperature dependent. More informa- 
tion on the neutrino spectrum, neutron-star mass and interior physics of the 
PNS will undoubtedly shed further light on the problem. However, the like- 
lihood function decreases rapidly as 02 and are increased (see Fig. 2 for 
0-2 > 0.01), and so the bounds on these quantities may not be greatly affected 
by advances in our knowledge of the physics which governs the propagation 
of energy in the PNS interior. 

In contrast, the bound derived here is sensitive to the rate at which energy 
is radiated to the GODs, i.e. to the temperature-dependence of the graviton 
emissivity. In particular, many-body effects, such as the multiple-scattering 
suppression of bremsstrahlung reactions which occurs in a dense plasma and 
is known as the Landau- Pomeranchuk-Migdal effect, might strongly reduce the 
KK-gravistrahlung energy-loss rate. Indeed, earlier studies [20] of the multiple- 
scattering suppression of axion radiation due to strong nucleon-spin fluctua- 
tions showed that the axion-mass bound was weakened by a factor of two 
when such effects were considered. Many-body effects might reasonably be 
expected to have a similar influence here — although gravitons couple to the 
stress-energy tensor, and not, as axions do, to the nucleon spin, and so they 
are sensitive to different correlations in the dense plasma. All of these issues 
warrant further investigation and they will be addressed in future work. In this 
work we have derived the bounds (8) and (9), which improve upon previous 
SN 1987a bounds on the size and scales of GODs in three important ways: 

(1) They represent the flrst self-consistent simulations of the early, neutrino- 
emitting phase of a proto-neutron star which also include energy losses 
due to the couphng of the Kaluza-Klein modes of the graviton have been 
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employed to examine the SN 1987a neutrino signal in detail. 

(2) They use the KK-graviton emissivities of Ref. [4], thereby anchoring the 
nuclear physics of the emissivity for the A^A^ KK-gravistrahlung on solid 
ground. 

(3) They are derived using a rigorous probabilistic analysis that facilitates an 
unambiguous criterion for the likelihood of certain values of the GODs' 
radii, and a consequent estimation of confidence levels. 

While these improvements do not alter the earlier estimated bounds signifi- 
cantly they do serve to place them on a firm theoretical basis. We also believe 
that in the event of a future galactic supernova the formalism and ideas devel- 
oped in this work will prove useful in analyses of the large numbers of counts 
expected in current-generation neutrino detectors, such as Super Kamiokande 
and SNO. 
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